Nuclear morphology is shaped by loop-extrusion programs.

It is well established that neutrophils adopt malleable polymorphonuclear shapes to migrate through narrow interstitial tissue spaces1-3. However, how polymorphonuclear structures are assembled remains unknown4. Here we show that in neutrophil progenitors, halting loop extrusion-a motor-powered process that generates DNA loops by pulling in chromatin5-leads to the assembly of polymorphonuclear genomes. Specifically, we found that in mononuclear neutrophil progenitors, acute depletion of the loop-extrusion loading factor nipped-B-like protein (NIPBL) induced the assembly of horseshoe, banded, ringed and hypersegmented nuclear structures and led to a reduction in nuclear volume, mirroring what is observed during the differentiation of neutrophils. Depletion of NIPBL also induced cell-cycle arrest, activated a neutrophil-specific gene program and conditioned a loss of interactions across topologically associating domains to generate a chromatin architecture that resembled that of differentiated neutrophils. Removing NIPBL resulted in enrichment for mega-loops and interchromosomal hubs that contain genes associated with neutrophil-specific enhancer repertoires and an inflammatory gene program. On the basis of these observations, we propose that in neutrophil progenitors, loop-extrusion programs produce lineage-specific chromatin architectures that permit the packing of chromosomes into geometrically confined lobular structures. Our data also provide a blueprint for the assembly of polymorphonuclear structures, and point to the possibility of engineering de novo nuclear shapes to facilitate the migration of effector cells in densely populated tumorigenic environments.

Single cell transcriptomic profiling of a neuron-astrocyte assembloid tauopathy model.

The use of iPSC derived brain organoid models to study neurodegenerative disease has been hampered by a lack of systems that accurately and expeditiously recapitulate pathogenesis in the context of neuron-glial interactions. Here we report development of a system, termed AstTau, which propagates toxic human tau oligomers in iPSC derived neuron-astrocyte assembloids. The AstTau system develops much of the neuronal and astrocytic pathology observed in tauopathies including misfolded, phosphorylated, oligomeric, and fibrillar tau, strong neurodegeneration, and reactive astrogliosis. Single cell transcriptomic profiling combined with immunochemistry characterizes a model system that can more closely recapitulate late-stage changes in adult neurodegeneration. The transcriptomic studies demonstrate striking changes in neuroinflammatory and heat shock protein (HSP) chaperone systems in the disease process. Treatment with the HSP90 inhibitor PU-H71 is used to address the putative dysfunctional HSP chaperone system and produces a strong reduction of pathology and neurodegeneration, highlighting the potential of AstTau as a rapid and reproducible tool for drug discovery.

FOXR1 regulates stress response pathways and is necessary for proper brain development.

The forkhead box (Fox) family of transcription factors are highly conserved and play essential roles in a wide range of cellular and developmental processes. We report an individual with severe neurological symptoms including postnatal microcephaly, progressive brain atrophy and global developmental delay associated with a de novo missense variant (M280L) in the FOXR1 gene. At the protein level, M280L impaired FOXR1 expression and induced a nuclear aggregate phenotype due to protein misfolding and proteolysis. RNAseq and pathway analysis showed that FOXR1 acts as a transcriptional activator and repressor with central roles in heat shock response, chaperone cofactor-dependent protein refolding and cellular response to stress pathways. Indeed, FOXR1 expression is increased in response to cellular stress, a process in which it directly controls HSPA6, HSPA1A and DHRS2 transcripts. The M280L mutant compromises FOXR1's ability to respond to stress, in part due to impaired regulation of downstream target genes that are involved in the stress response pathway. Quantitative PCR of mouse embryo tissues show Foxr1 expression in the embryonic brain. Using CRISPR/Cas9 gene editing, we found that deletion of mouse Foxr1 leads to a severe survival deficit while surviving newborn Foxr1 knockout mice have reduced body weight. Further examination of newborn Foxr1 knockout brains revealed a decrease in cortical thickness and enlarged ventricles compared to littermate wild-type mice, suggesting that loss of Foxr1 leads to atypical brain development. Combined, these results suggest FOXR1 plays a role in cellular stress response pathways and is necessary for normal brain development.

Single cell transcriptomics reveals opioid usage evokes widespread suppression of antiviral gene program.

Chronic opioid usage not only causes addiction behavior through the central nervous system, but also modulates the peripheral immune system. However, how opioid impacts the immune system is still barely characterized systematically. In order to understand the immune modulatory effect of opioids in an unbiased way, here we perform single-cell RNA sequencing (scRNA-seq) of peripheral blood mononuclear cells from opioid-dependent individuals and controls to show that chronic opioid usage evokes widespread suppression of antiviral gene program in naive monocytes, as well as in multiple immune cell types upon stimulation with the pathogen component lipopolysaccharide. Furthermore, scRNA-seq reveals the same phenomenon after a short in vitro morphine treatment. These findings indicate that both acute and chronic opioid exposure may be harmful to our immune system by suppressing the antiviral gene program. Our results suggest that further characterization of the immune modulatory effects of opioid is critical to ensure the safety of clinical opioids.

ATAC-seq Assay with Low Mitochondrial DNA Contamination from Primary Human CD4+ T Lymphocytes.

ATAC-seq has become a widely used methodology in the study of epigenetics due to its rapid and simple approach to mapping genome-wide accessible chromatin. In this paper we present an improved ATAC-seq protocol that reduces contaminating mitochondrial DNA reads. While previous ATAC-seq protocols have struggled with an average of 50% contaminating mitochondrial DNA reads, the optimized lysis buffer introduced in this protocol reduces mitochondrial DNA contamination to an average of 3%. This improved ATAC-seq protocol allows for a near 50% reduction in the sequencing cost. We demonstrate how these high-quality ATAC-seq libraries can be prepared from activated CD4+ lymphocytes, providing step-by-step instructions from CD4+ lymphocyte isolation from whole blood through data analysis. This improved ATAC-seq protocol has been validated in a wide range of cell types and will be of immediate use to researchers studying chromatin accessibility.

Genetic determinants of co-accessible chromatin regions in activated T cells across humans.

Over 90% of genetic variants associated with complex human traits map to non-coding regions, but little is understood about how they modulate gene regulation in health and disease. One possible mechanism is that genetic variants affect the activity of one or more cis-regulatory elements leading to gene expression variation in specific cell types. To identify such cases, we analyzed ATAC-seq and RNA-seq profiles from stimulated primary CD4+ T cells in up to 105 healthy donors. We found that regions of accessible chromatin (ATAC-peaks) are co-accessible at kilobase and megabase resolution, consistent with the three-dimensional chromatin organization measured by in situ Hi-C in T cells. Fifteen percent of genetic variants located within ATAC-peaks affected the accessibility of the corresponding peak (local-ATAC-QTLs). Local-ATAC-QTLs have the largest effects on co-accessible peaks, are associated with gene expression and are enriched for autoimmune disease variants. Our results provide insights into how natural genetic variants modulate cis-regulatory elements, in isolation or in concert, to influence gene expression.

Genome-scale identification of transcription factors that mediate an inflammatory network during breast cellular transformation.

Transient activation of Src oncoprotein in non-transformed, breast epithelial cells can initiate an epigenetic switch to the stably transformed state via a positive feedback loop that involves the inflammatory transcription factors STAT3 and NF-κB. Here, we develop an experimental and computational pipeline that includes 1) a Bayesian network model (AccessTF) that accurately predicts protein-bound DNA sequence motifs based on chromatin accessibility, and 2) a scoring system (TFScore) that rank-orders transcription factors as candidates for being important for a biological process. Genetic experiments validate TFScore and suggest that more than 40 transcription factors contribute to the oncogenic state in this model. Interestingly, individual depletion of several of these factors results in similar transcriptional profiles, indicating that a complex and interconnected transcriptional network promotes a stable oncogenic state. The combined experimental and computational pipeline represents a general approach to comprehensively identify transcriptional regulators important for a biological process.

PHF6 regulates phenotypic plasticity through chromatin organization within lineage-specific genes.

Developmental and lineage plasticity have been observed in numerous malignancies and have been correlated with tumor progression and drug resistance. However, little is known about the molecular mechanisms that enable such plasticity to occur. Here, we describe the function of the plant homeodomain finger protein 6 (PHF6) in leukemia and define its role in regulating chromatin accessibility to lineage-specific transcription factors. We show that loss of Phf6 in B-cell leukemia results in systematic changes in gene expression via alteration of the chromatin landscape at the transcriptional start sites of B-cell- and T-cell-specific factors. Additionally, Phf6KO cells show significant down-regulation of genes involved in the development and function of normal B cells, show up-regulation of genes involved in T-cell signaling, and give rise to mixed-lineage lymphoma in vivo. Engagement of divergent transcriptional programs results in phenotypic plasticity that leads to altered disease presentation in vivo, tolerance of aberrant oncogenic signaling, and differential sensitivity to frontline and targeted therapies. These findings suggest that active maintenance of a precise chromatin landscape is essential for sustaining proper leukemia cell identity and that loss of a single factor (PHF6) can cause focal changes in chromatin accessibility and nucleosome positioning that render cells susceptible to lineage transition.

Systematic Epigenomic Analysis Reveals Chromatin States Associated with Melanoma Progression.

The extent and nature of epigenomic changes associated with melanoma progression is poorly understood. Through systematic epigenomic profiling of 35 epigenetic modifications and transcriptomic analysis, we define chromatin state changes associated with melanomagenesis by using a cell phenotypic model of non-tumorigenic and tumorigenic states. Computation of specific chromatin state transitions showed loss of histone acetylations and H3K4me2/3 on regulatory regions proximal to specific cancer-regulatory genes in important melanoma-driving cell signaling pathways. Importantly, such acetylation changes were also observed between benign nevi and malignant melanoma human tissues. Intriguingly, only a small fraction of chromatin state transitions correlated with expected changes in gene expression patterns. Restoration of acetylation levels on deacetylated loci by histone deacetylase (HDAC) inhibitors selectively blocked excessive proliferation in tumorigenic cells and human melanoma cells, suggesting functional roles of observed chromatin state transitions in driving hyperproliferative phenotype. Through these results, we define functionally relevant chromatin states associated with melanoma progression.

Iterative Modeling Reveals Evidence of Sequential Transcriptional Control Mechanisms.

Combinatorial control of gene expression is presumed to be mediated by molecular interactions between coincident transcription factors (TFs). While information on the genome-wide locations of TFs is available, the genes they regulate and whether they function combinatorially often remain open questions. Here, we developed a mechanistic, rather than statistical, modeling approach to elucidate TF control logic from gene expression data. Applying this approach to hundreds of genes in 85 datasets measuring the transcriptional responses of murine fibroblasts and macrophages to cytokines and pathogens, we found that stimulus-responsive TFs generally function sequentially in logical OR gates or singly. Logical AND gates were found between NF-κB-responsive mRNA synthesis and MAPKp38-responsive control of mRNA half-life, but not between temporally coincident TFs. Our analyses identified the functional target genes of each of the pathogen-responsive TFs and prompt a revision of the conceptual underpinnings of combinatorial control of gene expression to include sequentially acting molecular mechanisms that govern mRNA synthesis and decay.

Semiconductor-based DNA sequencing of histone modification states.

The recent development of a semiconductor-based, non-optical DNA sequencing technology promises scalable, low-cost and rapid sequence data production. The technology has previously been applied mainly to genomic sequencing and targeted re-sequencing. Here we demonstrate the utility of Ion Torrent semiconductor-based sequencing for sensitive, efficient and rapid chromatin immunoprecipitation followed by sequencing (ChIP-seq) through the application of sample preparation methods that are optimized for ChIP-seq on the Ion Torrent platform. We leverage this method for epigenetic profiling of tumour tissues.

Mechanisms establishing TLR4-responsive activation states of inflammatory response genes.

Precise control of the innate immune response is required for resistance to microbial infections and maintenance of normal tissue homeostasis. Because this response involves coordinate regulation of hundreds of genes, it provides a powerful biological system to elucidate the molecular strategies that underlie signal- and time-dependent transitions of gene expression. Comprehensive genome-wide analysis of the epigenetic and transcription status of the TLR4-induced transcriptional program in macrophages suggests that Toll-like receptor 4 (TLR4)-dependent activation of nearly all immediate/early- (I/E) and late-response genes results from a sequential process in which signal-independent factors initially establish basal levels of gene expression that are then amplified by signal-dependent transcription factors. Promoters of I/E genes are distinguished from those of late genes by encoding a distinct set of signal-dependent transcription factor elements, including TATA boxes, which lead to preferential binding of TBP and basal enrichment for RNA polymerase II immediately downstream of transcriptional start sites. Global nuclear run-on (GRO) sequencing and total RNA sequencing further indicates that TLR4 signaling markedly increases the overall rates of both transcriptional initiation and the efficiency of transcriptional elongation of nearly all I/E genes, while RNA splicing is largely unaffected. Collectively, these findings reveal broadly utilized mechanisms underlying temporally distinct patterns of TLR4-dependent gene activation required for homeostasis and effective immune responses.

The specificity of innate immune responses is enforced by repression of interferon response elements by NF-κB p50.

The specific binding of transcription factors to cognate sequence elements is thought to be critical for the generation of specific gene expression programs. Members of the nuclear factor κB (NF-κB) and interferon (IFN) regulatory factor (IRF) transcription factor families bind to the κB site and the IFN response element (IRE), respectively, of target genes, and they are activated in macrophages after exposure to pathogens. However, how these factors produce pathogen-specific inflammatory and immune responses remains poorly understood. Combining top-down and bottom-up systems biology approaches, we have identified the NF-κB p50 homodimer as a regulator of IRF responses. Unbiased genome-wide expression and biochemical and structural analyses revealed that the p50 homodimer repressed a subset of IFN-inducible genes through a previously uncharacterized subclass of guanine-rich IRE (G-IRE) sequences. Mathematical modeling predicted that the p50 homodimer might enforce the stimulus specificity of composite promoters. Indeed, the production of the antiviral regulator IFN-ß was rendered stimulus-specific by the binding of the p50 homodimer to the G-IRE-containing IFNß enhancer to suppress cytotoxic IFN signaling. Specifically, a deficiency in p50 resulted in the inappropriate production of IFN-ß in response to bacterial DNA sensed by Toll-like receptor 9. This role for the NF-κB p50 homodimer in enforcing the specificity of the cellular response to pathogens by binding to a subset of IRE sequences alters our understanding of how the NF-κB and IRF signaling systems cooperate to regulate antimicrobial immunity.

A unifying model for the selective regulation of inducible transcription by CpG islands and nucleosome remodeling.

We describe a broad mechanistic framework for the transcriptional induction of mammalian primary response genes by Toll-like receptors and other stimuli. One major class of primary response genes is characterized by CpG-island promoters, which facilitate promiscuous induction from constitutively active chromatin without a requirement for SWI/SNF nucleosome remodeling complexes. The low nucleosome occupancy at promoters in this class can be attributed to the assembly of CpG islands into unstable nucleosomes, which may lead to SWI/SNF independence. Another major class consists of non-CpG-island promoters that assemble into stable nucleosomes, resulting in SWI/SNF dependence and a requirement for transcription factors that promote selective nucleosome remodeling. Some stimuli, including serum and tumor necrosis factor-alpha, exhibit a strong bias toward activation of SWI/SNF-independent CpG-island genes. In contrast, interferon-beta is strongly biased toward SWI/SNF-dependent non-CpG-island genes. By activating a diverse set of transcription factors, Toll-like receptors induce both classes and others for an optimal response to microbial pathogens.

Epigenetic control: slow and global, nimble and local.

The regulation of gene expression involves multiple levels of control, from those that are inheritable to those that are highly responsive to environmental changes. In this issue of Genes & Development, Dong and colleagues (pp. 1159-1173) demonstrate that the dynamically controlled immune response transcription factor NF-kappaB may, in fact, have a role in regulating heterochromatin and gene expression at large distances from its actual target sequences and genes.

Sensitive ChIP-DSL technology reveals an extensive estrogen receptor alpha-binding program on human gene promoters.

ChIP coupled with microarray provides a powerful tool to determine in vivo binding profiling of transcription factors to deduce regulatory circuitries in mammalian cells. Aiming at improving the specificity and sensitivity of such analysis, we developed a new technology called ChIP-DSL using the DNA selection and ligation (DSL) strategy, permitting robust analysis with much reduced materials compared with standard procedures. We profiled general and sequence-specific DNA binding transcription factors using a full human genome promoter array based on the ChIP-DSL technology, revealing an unprecedented number of the estrogen receptor (ERalpha) target genes in MCF-7 cells. Coupled with gene expression profiling, we found that only a fraction of these direct ERalpha target genes were highly responsive to estrogen and that the expression of those ERalpha-bound, estrogen-inducible genes was associated with breast cancer progression in humans. This study demonstrates the power of the ChIP-DSL technology in revealing regulatory gene expression programs that have been previously invisible in the human genome.

Histone methylation-dependent mechanisms impose ligand dependency for gene activation by nuclear receptors.

Nuclear receptors undergo ligand-dependent conformational changes that are required for corepressor-coactivator exchange, but whether there is an actual requirement for specific epigenetic landmarks to impose ligand dependency for gene activation remains unknown. Here we report an unexpected and general strategy that is based on the requirement for specific cohorts of inhibitory histone methyltransferases (HMTs) to impose gene-specific gatekeeper functions that prevent unliganded nuclear receptors and other classes of regulated transcription factors from binding to their target gene promoters and causing constitutive gene activation in the absence of stimulating signals. This strategy, based at least in part on an HMT-dependent inhibitory histone code, imposes a requirement for specific histone demethylases, including LSD1, to permit ligand- and signal-dependent activation of regulated gene expression. These events link an inhibitory methylation component of the histone code to a broadly used strategy that circumvents pathological constitutive gene induction by physiologically regulated transcription factors.